Vortex— vortex interaction in two-component Bose— Einstein condensates 
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The vortex-vortex interaction in two-component Bose-Einstein condensates is shown to present 
characteristic effects not possible in single-component condensates. In particular, vortices in differ- 
ent components undergo separate, but concentric, orbits for both same and opposite circulation. In 
addition, contrary to the well-known behavior in single-component condensates, vortices with the 
same circulation created in the less dense component merge in the presence of dissipation, allowing 
for the creation of metastable vortices with more than one quantum of circulation. 



The experimental achievement of Bose-Einstein con- 
densation (BEC) in trapped dilute alkali gases ^ has 
opened the possibility to explore fundamental properties 
of low temperature physics in a very controllable way. In 
particular, superfluid phenomena, well known from ^He 
physics 1^, have been experimentally observed in BEC 
1^. Among the phenomena related with superfluidity, 
the possibility of vortices with quantized circulation 
has aroused a big interest. Due to the diluteness of the 
atomic condensates, accurate theoretical predictions are 
possible. In addition, in these systems, the diameter of 
the vortex core is typically three orders of magnitude 
larger than for superfluid Helium, allowing for a more 
straightforward optical observation. Numerous theoreti- 
cal schemes have been suggested for generating vortices 
in BEC, see e.g. Refs. P-pl|. So far only two tech- 
niques have successfully generated vortices in conden- 
sates, namely the stirring with a blue detuned laser [ p^ , 
and the coherent interconversion between the two com- 
ponents of a binary BEC |l3| . In addition the creation of 
a vortex ring after the decay of a dark soliton has been 
recently reported |Q . The dynamics and stability of vor- 
tices in trapped Bose-Einstein condensates has recently 
also been a subject of active investigation |p^ . 

In the recent years, the development of trapping tech- 
niques has allowed for the creation of multi-component 
condensates, which are formed by trapping atoms in dif- 
ferent internal (electronic) states The multicom- 
ponent BEC, far from being a trivial extension of the 
single-component one, presents novel and fundamentally 
different scenarios for its ground state 1 18 1^ and exci- 
tations 1^^. In particular, it has been observed that the 
condensate can reach an equilibrium state characterized 
by the separation of the species in different domains . 
The realization of vortices in binary condensates in re- 
cent experiments at JILA has aroused the interest for 
the properties of vortices in multicomponent condensates 
| p5| , pT|j22| . Several interesting phenomena have been re- 
ported in this context, such as vortex ground states ||2|] 
and generation of persistent currents ||2^] . 

In this Letter we discuss the interaction of vortices in 



a two-component BEC. The first part of the Letter is de- 
voted to homogeneous miscible binary condensates. We 
consider in particular the case in which one vortex is cre- 
ated in each component. In single-component BEC (or 
classical gas) in absence of dissipation, two vortices with 
the same circulation orbit around each other, whereas 
vortices with opposite circulation move parallelly (vortex 
pair) However, for the case of asymmetric nonlin- 
ear interactions between the two components of a binary 
BEC, the velocity fields for the vortices are asymmetric. 
We show that in this scenario, vortices with the same 
circulation perform two concentric orbits. Interestingly, 
also vortices with opposite circulation perform concentric 
orbits, clearly differing from the one component behavior. 

In the second part of this Letter we analyze the inter- 
action of two vortices created in one of the components of 
an inhomogeneous binary BEC. In particular, we study 
the interaction of vortices with the same circulation, and 
show that the "buoyancy" effect strongly distorts 
the vortex dynamics. In single-component BEC, vortices 
with the same circulation repel each other in the presence 
of dissipation, whereas the opposite can be true in a two- 
component BEC. This result strongly suggests that the 
creation of stable vortices with more than one quantum 
of circulation should be possible in a binary BEC. 

Let us first consider the case of one vortex created in 
each of the components of a miscible binary BEC trapped 
in a 2D disc-shape box-like potential of radius R. The 
BEC can be considered homogeneous, except in a re- 
gion of the size of the healing length, close to the box 
boundaries. The dynamics can be described in 2D if the 
mean-field interaction is smaller than the typical energy 
separation E±_ in the third direction To generate 

the vortices we use the phase-imprinting (PI) technique 
although other methods could be employed, as for 
example stirring methods [l^ , or the method employed 
in Ref. |[T^], which already involves two components. 

For sufficiently low temperatures the dynamics is well 
described by two coupled Gross-Pitaevskii equations 



«^^V'i(r,i) = 



1 



gji\ipi{r,t)\'^} iljj{r,t), 



(1) 



where V{r) is the trap potential, gji — 'inh^aji/mL is the 
coupling constant (averaged over the frozen-out transver- 
sal direction) between the components j and / with the 
transversal length L, m is the atomic mass and Uji the 
scattering length between j and I = 1,2). 

We have numerically simulated the creation of one vor- 
tex in each component by using PI. The creation of a vor- 
tex in one component modifies the density in the other 
one, due to the coupling in Eq. (^. For example, dark 
solitons in homogeneous two-component BEC can only 
be created at distances comparable to the healing length 
[p6t , Iq = h/y/gnm, where n is the density of the sam- 
ple. The latter does not apply to vortices, which can be 
created at distances significantly larger than Iq, without 
being destroyed by the density fluctuations of the other 
component. Although we have numerically observed sta- 
ble creation of vortices at distances larger than 20Zo, typi- 
cal distances of the order of 81q are used in the simulations 
below, since the time scale of the vortex-vortex dynamics 
is given by m(P /h, where d is the distance between the 
vortices. Therefore shorter distances allow for faster dy- 
namics. This facilitates the numerical simulations, and 
reduces instabilities during the vortex interaction. As al- 
ready observed in recent experiments , the vortex core 
in one component is filled by the other component. Al- 
though the two components interact through the density, 
and therefore they are not directly experiencing the cir- 
culation of the other vortex, the filling of each core does 
feel the circulation created by the vortex in the same 
component, dragging in its way the vortex core in which 
it is placed. The result is an indirect interaction. 

Let us postpone the discussion about the boundary ef- 
fects and concentrate on the situation in which the radius 
R of the disc-shape trap is sufficiently large to neglect the 
boundary effects during the time scale of the simulations. 
In the following we consider gn = g22 = 1.055i2. These 
particular g^j values have been chosen for simplicity, to 
guarantee a complete overlap between the wavefunctions 
of both components. In the case of *^Rb, the coupling 
constants have the values gn : 522 '■ 512 = 0.97 : 1.03 : 
1.0. These values also imply a large overlap, without 
separate domains, and therefore the dynamics we discuss 
below should be observable in *^Rb. 

In the symmetric situation, ni_2 = no, the qualitative 
behavior of single-component vortex-vortex interaction 
is recovered, namely, vortices with the same circulation 
rotate around each other, forming a vortex orbit, whereas 
vortices of opposite circulation move paralelly (vortex 
pair) through the fluid. The time evolution of the vortex 
cores is, however, different to the single-component one, 
since the effective masses of the vortex cores are different. 

This picture breaks down when the ni 7^ ^2. Fig. 
|l| shows the position of the vortex lines for the case of 
ni = 2n2 = lO^^cm"^, d = 4/im, and R = 20d, for 



vortices with the same circulation. In this asymmetric 
scenario, contrary to the case of single-component vor- 
tices, vortices with the same circulation do not undergo 
an orbit, but two concentric ones. The inner orbit corre- 
sponds to the vortex in the denser component, whereas 
the outer one corresponds to the vortex created in the less 
dense one. Interestingly, such a behavior is also present 
for the case of vortices with opposite circulation. 

In order to understand the formation of the planetary 
orbits when rii 7^ n2, let us consider the situation of a 
vortex line in the component 1 (2) at a position ri (f^). 
The density of the component 1 (2) is considered constant 
and equal to ni (712) everywhere in the system, except 
in the core of the vortices, which will be modeled as a 
cylinder of radius h (Z2). In the core of the vortex in the 
component 1 (2) the density of 1 (2) is considered zero, 
whereas the density of 2 (1) is considered ni -\- n2, in 
such a way that the total density is flat. In the following 
calculation, we shall assume that \r2\ ^ R, i.e. the 
vortices are placed close to the trap center, and the vortex 
cores are far away from each other, Zi , Z2 ^ 1^2 ^ | = H ■ 
The conservation of linear momentum and energy in the 
system, up to order 0((rj/i?)2), 0{il,/r)^), (j = 1,2), 
lead respectively to the following equations: 
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where the velocity of the vortex in the j-th component is 



given by Vj = uj / r 



U2 — U1, and K = 1 (—1) for vor- 



tices with the same (opposite) circulation. Neglecting in 
Eq. (||) the relatively small rhs, and assuming R tending 
to 00, one obtains the solution V2 ~ (ri — r2)"'"/|ri — r2p, 
and vi — — {Kn2 / ni)v2 l27| . The motion of the vor- 
tex lines is then determined by the first-order equations 



dfj I dt 



The trajectories provided by these equa- 



tions (solid lines in Fig. gj) describe very well the nu- 
merically calculated vortex orbits. Note that since the 
velocity field is always tangential to the vector which 
joins the vortex lines, the distance between the vortices 
remains constant at any time, and consequently the an- 
gular frequency for both orbits is the same. Therefore, 
vortices of the same circulation are always at opposite an- 
gles of their respective orbits, whereas vortices with op- 
posite circulation always occupy the same angle. Hence, 
the angular frequency for both orbits is the same. From 
this fact, it is easy to show that for vortices with the 
same/opposite circulation, initially placed at a distance 
d, the orbit of the outer vortex (say in component 1) has 
a radius r\ = dn2/{ni ±712), whereas the inner one has 
a radius r2 — dni/irii ±712). Note that for rii = 712 the 
familiar one-component results are retrieved. 

Let us discuss at this point the role of the boundaries 
on the dynamics. Our numerical analysis shows that the 
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finite size of the trap does not disturb the distance be- 
tween the vortices, and, as in the case of i? ^ oo, r i — r2 
moves in a circle. However, the coordinate (ri + '^2)/2, 
which for R oo moves in a circle as well, performs 
in the presence of a finite R a cycloid motion. As a re- 
sult, the planetary orbits are not circular anymore, but 
will precess around a circular orbit, as one can observe 
in Fig. H, where the same parameters of Fig. |l| are con- 
sidered, except R = 4.5d. Nevertheless, similar to the 
case R ^ oo, the trajectory of the vortex in the denser 
component is always the inner trajectory. 

In the last part of this Letter we analyze the case of a 
binary BEC trapped in a 2D harmonic potential of fre- 
quency Lo. This resembles the situation in which vortices 
are currently created in *^Rb in JILA In particu- 

lar, we shall analyze the interaction of two vortices in the 
same component, which are created symmetrically from 
the condensate center, and with the same circulation. In 
homogeneous as well as single-component trapped con- 
densates, and in the presence of dissipation, the vortices 
will spiral out from their initial orbit until they finally 
reach the boundaries of the condensate cloud, where they 
are destroyed . We shall show that this is not the case 
in an inhomogeneous binary BEC. 

In a recent paper ||2^, the dynamics of a filled vor- 
tex in a trapped binary BEC was analyzed. In an 
inhomogeneous potential, the lowest energy solution is 
reached when the component with smaller nonlinear self- 
interaction occupies the regions of highest density, i.e. 
the trap center. Therefore, in our case, if the less dense 
component is displaced from the trap center, it will tend 
to return to it. This effect is sometimes called the "buoy- 
ancy" effect This effect can counteract the ten- 
dency of an off-centered vortex to spiral out. In par- 
ticular, it has been shown that for sufficiently large non- 
linearity a vortex, created in the component with larger 
self-interaction at the center of the trap, can become 
metastable |2^. We employ a similar treatment to ana- 
lyze the interaction of vortices. It consists in evaluating, 
for different distances between vortices, the free energy 

I d3r(|l(|V^i(r)|2 + |VV^2(r)|2) 

+Vi{r)\Mr)\^ + V2{r)\Mr)\^ 

+ ^\Mrt+9i2\Mr)\'\Mr)\' + ^\Mr)\'')- 

In order to obtain a clear value of the free energy associ- 
ated with every vortex configuration, instead of creating 
the vortices using the PI method discussed above, we 
generate them, as in Ref. [Q, by evolving the Gross- 
Pitaevskii equation (|l|) in imaginary time, fixing at ev- 
ery iteration the phase of the component in which the 
vortices are created, in order to have the correct phase 
properties. Once imprinted in this way two vortices at 
position (— fi/2,0) and ((i/2,0), we calculate the energy 



E{d) following Eq. (||). For the case of single component 
condensates, one can show that E{d) decreases with d, 
for any d. This reflects the fact that in the presence of 
dissipation the vortices tend to spiral out. For vortices 
created in the denser component (say 1), the situation is, 
however, completely different. Fig. ^ shows the case of 
^^Rb (ail — 1.03ai2, 022 — 0.97ai2, and 012 = 5.7nm), 
with a trap frequency uj — 2tt x SOHz and L ~ 1/im, with 
number of particles iVi — 10^ and iV2 = 8 x 10^. One 
observes that E{d) clearly grows with d for low d, until it 
reaches a maximal value. This maximum corresponds to 
the case in which the vortices are placed at the maxima of 
the unperturbed density distribution of the component 1 
(see inset in Fig. ||) . For larger distances the tendency of 
the vortices to spiral out is no more counteracted since 
the slope of the density distribution changes the sign. 
Therefore, vortices created sufficiently close to the center 
will tend to spiral inwards in the presence of dissipation, 
until eventually merging at the trap center. Once they 
merge, a vortex with double circulation will be created. 
Since due to its topology a vortex can only be destroyed 
if it reaches a region of zero density, the created vor- 
tex of double circulation will remain metastable at the 
trap center, where only inelastic processes will in prac- 
tice limit its lifetime. Following similar arguments, the 
merging of vortices should allow for the creation of a 
metastable vortex with n.^ quanta of circulation, open- 
ing a new scenario which is naturally forbidden in both 
homogeneous single- and multi-component BEC, as well 
as in inhomogeneous single-component condensates. 

Let us finally point out, that a similar argument ap- 
plies for vortices with opposite circulation created in the 
denser component and moving parallelly in a toroidal 
trap In that case, the "buoyancy" effect counteracts 
the natural tendency of the vortices to mutually annihi- 
late, allowing for the creation of a vortex pair, which will 
circulate around the torus center in a metastable orbit. 

Summarizing, we have shown that the vortex-vortex 
interaction in multicomponcnt condensates, presents 
characteristic features such as planetary orbits for vor- 
tices placed in different components, and opens new 
interesting possibilities for creating metastable vortices 
with more than one quanta of circulation, as well as 
metastable vortex pairs in toroidal geometries. 

The analysis of the vortex-vortex interaction has been 
limited to the case of 2D traps. Lower dimensional con- 
densates are currently actively investigated, and very re- 
cently the experimental observation of 2D (and even ID) 
condensates have been reported |^,^. However, the 
results of the present Letter, open interesting questions 
concerning the behavior of vortices in 3D binary BEC. 
In particular, for linear vortices, if the relative density 
between the different components varies in the direc- 
tion of the vortex lines, it should be expected that vor- 
tices created in different components will undergo at each 
transversal position planetary orbits of different periods 
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and radii. Consequently, non trivial distortions of the 
vortex lines are expected. In addition, the recent obser- 
vation of vortex-rings in condensates stimulates the 
investigation of the interaction of such more complicated 
structures in binary BEC. The analysis of these effects is 
beyond the scope of the present Letter, and will be the 
subject of further investigation. 
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FIG. 1. Numerical (circles) and analytical (solid lines) re- 
sults for the vortex orbits with ni — 2n2, d = 4/im, R = 20d. 
The denser component occupies the inner orbit. 
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FIG. 2. Same as Fig. m but with R = 4.5d. 
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FIG. 3. Energy as a function of vortex distance (see text 
for details). Inset: Density distribution without vortices. 
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